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ABSTRACT: A model was developed to assess the impact of chromophoric dissolved organic matter 
(CDOM) on phytoplankton production within the euphotic zone. The rate of depth-integrated daily 
gross primary productivity within the euphotic zone (^GPR.*) was evaluated as a function of date, 
latitude, CDOM absorption (a< DOM ) characteristics, chlorophyll a (chi a) concentration, vertical stratifi- 
cation, and phytoplankton sensitivity to UV radiation (UVR). Results demonstrated that primary pro- 
duction was enhanced in the upper -30 m of the water column by the presence of CDOM, where pre- 
dicted increases in production due to the removal of damaging UVR more than offset its reduction 
resulting from the absorption of photosynthetically usable radiation. At greater depths, where little 
UVR remained, primary production was always reduced due to removal by CDOM of photosyntheti- 
cally usable radiation, When CDOM was distributed homogeneously within the euphotic zone, J z GPP (l/ 
was reduced under most bio-optical (i.e. solar zenith angle, chi a and CDOM absorption, and ozone 
concentration) and photophysiological (i.e. sensitivity to UVR) conditions because the predicted reduc- 
tion in primary production at depth was greater than the enhancement of production at the surface. 
A reduction in j^GPP,,, was also predicted when CDOM or phytoplankton was restricted to near-surface 
waters (-30 m) and CDOM absorption was moderate [a ( [)nM (450) > 0.015 m ']. j„GPP ( , /( however, was 
enhanced when CDOM or phytoplankton was restricted to a very shallow surface layer (-10 m), even 
if CDOM absorption was high |<3 (T)om(^) at 450 nm -0.07 m *]. Changes in J^GPP^ resulting from the 
presence of CDOM were only slightly sensitive to ozone concentrations. In well-mixed waters where 
the flux of UVB is relatively high, such as in the Southern Ocean when the ozone hole is present, the 
presence of CDOM should result in little or no enhancement of j^GPP,,,, although phytoplankton pro- 
duction would be expected to increase somewhat in surface waters. 
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INTRODUCTION 

Recent studies have demonstrated that ambient 
levels of ultraviolet radiation (UVR, 280 to 400 nm) can 
substantially decrease rates of carbon fixation by 
phytoplankton (Karentz et al. 1991, Cullen et al. 1992, 
Helbling et al. 1992, Smith et al. 1992, Holm-Hansen 
et al, 1993a). At normal ozone (0 3 ) concentrations, i.e. 
344 Dobson Units (DU), UVR has been found to reduce 
primary productivity in surface waters by as much as 
50% (Cullen et al. 1992, Holm-Hansen et al. 1993b, 
Cullen & Neale 1994). Because seawater strongly 
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attenuates UVR, particularly at shorter wavelengths, 
UVR inhibition is most severe near the sea surface. For 
example, along 50° S in mid December, a normal level 
of UVR was estimated to reduce phytoplankton pro- 
duction by 57 % at a depth of 1 m, while UVR inhibition 
decreased to <5% at 30 m (Arrigo 1994). Inhibition 
associated with UVR is due primarily to UVA (320 to 
400 nm), which is more abundant at the sea surface 
and is attenuated less strongly in the water column 
than the more damaging UVB (280 to 320 nm). Esti- 
mates of UVR inhibition resulting from 0 3 depletion 
range from an additional 1 to 12% of depth-integrated 
daily primary productivity (Smith et al. 1992, Holm- 
Hansen et al. 1993a, Arrigo 1994), due to the increased 
atmospheric transmission of UVB (UVA is not affected 
by changes in 0 3 ). 


208 


Mar Ecol Prog Ser 140: 207-216, 1996 


The degree to which the increased flux of incoming 
UVR affects phytoplankton photosynthesis depends in 
part upon the optical characteristics of the water col- 
umn. In addition to pure seawater, which attenuates 
radiation between 280 and 400 nm, particulates such 
as phytoplankton, bacteria, and detritus also absorb 
UVR. Holm-Hansen et al. (1993b) measured total par- 
ticulate absorption coefficients at 300 nm as ranging 
from 0.1 to 0.3 m 1 in Southern Ocean waters. Specific 
absorption coefficients for Antarctic phytoplankton 
have been reported to exceed 0.1 m 2 (mg chi a) -1 
within the UV range (Mitchell et al. 1989, Arrigo 1994). 
Kopelevich et al. (1987) showed that bacterial attenua- 
tion at 390 nm ranged from 0.002 m 1 for Micrococcus 
sp. to 2.80 m 1 for Moraxella sp. at concentrations of 
1012 cells nr 3 , and increased markedly at shorter 
wavelengths. 

Many dissolved constituents within seawater, such 
as chromophoric dissolved organic matter (CDOM), 
are also strong absorbers of UVR, often having absorp- 
tion coefficients that are higher than those of both sea- 
water and suspended particulates. CDOM absorption 
within the UV is greatest at 280 nm and declines expo- 
nentially at greater wavelengths (Fig. 1). The slope of 
log(absorption) versus wavelength is typically in the 
range of 0.012 to 0.029 nm' 1 (Bricaud et al. 1981, 
Carder et al. 1989, Blough et al. 1993, Hoge et al. 1993, 
Green & Blough 1994). 

CDOM absorption can be highly variable both tem- 
porally and spatially, changing as a function of CDOM 
concentration and specific absorption (Carder et al. 
1989, Blough et al. 1993). Carbon-specific absorption 
coefficients for riverine sources of CDOM are 10 to 150 
times greater than those of marine CDOM (Carder et 
al. 1989), either because these forms of CDOM are 
more strongly absorbing or because a greater propor- 
tion of riverine DOC consists of CDOM. Consequently, 
coastal waters typically possess higher CDOM absorp- 



Fig. 1 . Equation used to describe the carbon-specific absorption 
spectrum for chromophoric dissolved organic matter (^cdom) 
and its corresponding shape. The coefficient acDOM(450) = 
0.03 m 1 and the spectral slope S = 0.015 nm" 1 


tion coefficients than open ocean waters. For example, 
CDOM absorption in the near-shore waters of the 
Delaware Bay, USA, ranged from <0.05 to 0.38 m 1 at 
442 nm (M. D. DeGrandpre unpubl.) while those mea- 
sured at offshore sites in the Gulf of Mexico were 
approximately an order of magnitude lower, varying 
from 0.002 to 0.068 m 1 at 440 nm (Carder et al. 1989, 
Green & Blough 1994). 

As an absorber of radiation within the water column, 
CDOM plays 2 distinct roles in the marine environ- 
ment. First, because CDOM strongly absorbs UVR, it 
may act as a UV photoprotectant and reduce the dele- 
terious effect of UVR on phytoplankton, particularly 
in regions of ozone depletion, such as the Southern 
Ocean, where fluxes of UVR are increasing. Second, 
CDOM also attenuates photosynthetically usable radi- 
ation (PUR, that portion of photosynthetically available 
radiation, PAR, between 400 nm and 700 nm that can 
potentially be absorbed by phytoplankton pigments), 
which will reduce primary productivity in regions 
where light is limiting. The primary objective of this 
study was to assess the relative importance of these 2 
opposing influences of CDOM and evaluate how these 
processes vary under a wide variety of environmental 
conditions. Specifically, we investigated how the im- 
pact of CDOM on depth-integrated rates of primary 
production changes as a function of the vertical distri- 
butions and concentrations of chi a and CDOM, CDOM 
absorption properties, latitude, and stratospheric ozone 
concentration. We have restricted our study to the 
open ocean where CDOM absorption is less variable 
and therefore more easily characterized, though our 
general conclusions are likely to be applicable in 
coastal waters. 


METHODS 

The 1 -dimensional model of Arrigo (1994) was used 
to compute the depth-dependent variation in primary 
production as a function of the spectral distribution of 
incoming solar radiation (280 to 700 nm). Included in 
the model were separate components for calculating 
the flux of atmospheric radiation, in-water bio-optics, 
and primary production. The model calculated the rate 
of gross primary production (mg C m -3 d" 1 ) at 1 m 
vertical resolution on an hourly basis as a function of 
diurnal changes in spectral irradiance. UV inhibition of 
primary production was calculated using the biological 
weighting function (e) provided by Cullen et al. (1992). 
Production was integrated over depth (the euphotic 
depth, 0.1% of surface PAR) and time (24 h) to deter- 
mine daily production on an areal basis (J z GPP ez , mg C 
m~ 2 d -1 ). Note that J z GPP ez is defined here as the depth 
integrated rate of primary productivity, not the produc- 
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tivity at a given depth. Model input included day of the 
year, latitude, seawater temperature (°C), stratospheric 
[0 3 ] (DU), vertical distributions of phytoplankton (mg 
C m -3 ) and CDOM (g CDOM-C m“ 3 ) concentrations, 
and biomass-specific absorption coefficients. The 
effects of changes in UVR and PUR were determined 
independently by setting a flag that controls whether 
UV inhibition was calculated or not. If inhibition was 
not calculated, changes in production resulting from 
diminished 0 3 were due to variation in PUR only. 

Details of the model used in the present analysis can 
be found in Arrigo (1994); for convenience, the basic 
formulations are listed in Table 1. Diffuse attenuation 
by CDOM, not modeled by Arrigo (1994), has been 
described by the equation 


^cdom(^) 


a CDOM(^) 

P 


( 1 ) 


where K CDOM (X) a nd ^cdomM are the diffuse attenua- 
tion and absorption coefficients (nr 1 ), respectively, for 
CDOM, and p (dimensionless) is the mean cosine of the 
angular distribution of the downwelling light field. It 


has been shown that <3cdom(>M takes the form 

a cn omW - a cDOM(A. r )exp[S(X-X r )l (2) 

where a C DOM(^r) is the absorption by CDOM at the 
reference wavelength, ^ r , and S (nm' 1 ) is the slope of 
log(absorption) versus wavelength. 

A total of 512 simulations were performed. The con- 
ditions for the standard simulations were temperature 
= 0°C, X r = 450 nm, a CD o M (450) = 0.03 nr 1 , S = 0.015 
nm -1 , and chi a = 0.1 mg m” 3 . Ozone concentrations 
tested were 172 and 344 DU, latitude ranged from 
45° S to 75° S, and day of the year varied from 220 to 
360. UV inhibition of primary production was routinely 
included in the standard simulation. All quantities 
were assumed to be constant with depth unless other- 
wise specified. The values for 3 cdom( 450) and S used 
in Eq. (2) were derived by averaging similar measure- 
ments presented in Bricaud et al. (1981) and Green & 
Blough (1994) for the open ocean. In most cases, sensi- 
tivity analyses were compared to results of the stan- 
dard simulation and the differences between 2 corre- 
sponding runs are presented. 


Table 1. Equations used in the model. Details can be found in Arrigo (1994) and Arrigo & Sullivan (1994) 


Spectral light distribution 
(pEin m 2 s' 1 ): 

Bulk diffuse attenuation 
coefficient (m 1 )-. 

Diffuse attenuation for 
particulates (nr 1 ): 

Mean cosine for 
downwelling radiation: 

Gross microalgal production 
(mg C m“ 3 h -1 ): 

Rate of phytoplankton 
growth (d 1 ): 

Maximum phytoplankton 
growth rate (d _1 ): 

Light limitation 
(dimensionless): 


E (i {X,z) = E d (X,z-Az) exp[-A"(Vz)Az] 

K{X,z) = K W (X) -f K v (X,z) 

rs n WUl + a'cWQz) 

KplA,Z) = 

u = cos(0 w ) for direct light, 0,83 for diffuse 

<*(z,t) = G{z,t)P{z,t ) 

(H 

G{z,t) = p(z,t)G md]I (z) 


G m „ x (z) = G„ exp[rT(z)) 

p(z) = [~1 - expf - PUR(Z 1 gll f ) 

L h / fc '(z) V JJU + F‘, nh (z)J 


Photoadaptation parameter 
(pEin m 2 s' 1 ): 


j • i z \ - ^ max 

k 1 + 2exp|-B PUR'] 

B= exp[l. 089 -2.12 log (7 t ' mdX )] 


Mean PUR above D E during 
F (pEin m 2 s" 1 ): 


PUR’ = 


r\2+F/2 rl) v 

\ I PUR(z,f)dz df 

_ -Ui 2 r /2 

r12*F/2 f O E 


rV2*f/2 j *o E 

i=12-F/2 J/-0 


dzdt 


C:chl a at PUR(z,t) and I k ’(z): 

Photosynthetically usable 
radiation (pEin m' 2 s' 1 ): 

UV inhibition (dimensionless): 


_ 33,125 + 9.77F PUR(z) 
f ” 33. 125 + 9. 77F-/ k (z) 

PUR(z) = f 00 E d (X,z) a ' c dl 

' Jx-400 ' a* 

a cmax 

f400 

^‘inh(z) = L 28(1 E <^- E 'l(^.z)dX 


X (nm) is wavelength 
z (m) is depth 
t (h) is time 

K w (m 1 ) is the extinction coefficient for 
clear ocean water 

b ( (m 1 ) is the backscatter by micro- 
algae 

a’ r (m 2 mg chi a 1 ) is the pigment- 
specific absorption coefficient for 
phytoplankton 

C (mg chi a m 3 ) is the pigment con- 
centration 

P (mg C m H ) is the phytoplankton 
standing crop 

G () (d 1 ) is the microalgal growth rate at 
0°C 

T (°C) is temperature 
r (°C _1 ) is a rate constant that deter- 
mines the sensitivity of G mdX to T 

Vmdx is the maximum observed value 
for I k ' 

F (h) is the photoperiod 
D e (m) is the euphotic depth 

a* rmdX (m 2 m 9 chi a_1 ) is the maximum 
value attained by a\{X) 
e{X) [(mW m -2 )' 1 ] is the wavelength- 
dependent efficiency for damage to 
photosynthesis by UV 
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Day of the Year 

Fig. 2, Seasonal changes in gross primary productivity (mg C 
m“ 2 d 1 ) as a function of latitude for the standard simulations. 
Each data point represents a separate simulation 


RESULTS 

In the standard run, primary productivity varied as a 
function of the day of the year and latitude, ranging 
from 0 to 286 mg C m“ 2 d~ l (Fig. 2). Productivity was 
lowest in late austral winter at 75° S when the photo- 
period was <1 h and the sun remained low in the sky. 
Rates of production increased with time at all latitudes, 
the increase being most pronounced at high latitudes 
where changes in photoperiod were most rapid. Pre- 
dicted rates of production were highest during the 
summer solstice at 75° S, the time and location of max- 
imum photoperiod (24 h). Although the solar elevation 
was greater at lower latitudes, this did not compensate 
for the longer photoperiod at higher latitudes. 


To assess the amount of UV inhibition of j z GPP ez 
under a given suite of conditions (i.e. ozone concentra- 
tion, day of the year, and latitude), we calculated the 
difference between 2 corresponding simulations, one 
where UV inhibition was calculated (standard simula- 
tion) and the other where it was neglected. The spatial 
and temporal variation in UV inhibition of J z GPP ez 
was determined for a C DOM(450) = 0 and a C DOM(450) = 
0.03 m l as ozone concentration varied from 172 to 
344 DU (Fig. 3A, C). Results show that in the absence 
of CDOM, UV inhibition reduced J z GPP ez by as much 
as 9.2%. Inhibition was greatest at low latitudes in 
early austral summer (Fig. 3A). This was due to the 
increased flux of UVR resulting from the shorter path- 
length through the atmosphere at these increased solar 
elevations. When CDOM was present at all depths, 
model results followed the same pattern but UV inhibi- 
tion of j z GPP ez was reduced by -50% (Fig. 3B). The 
magnitude of this reduction is a measure of the mitiga- 
tion of UV inhibition provided by the presence of 
CDOM, i.e. its UV photoprotective capacity (Fig. 4A). 
The level of UV inhibition of J z GPP ez , and of photopro- 
tection by CDOM, changed only slightly when ozone 
concentration was reduced by 50% (Fig. 3C, D). 

CDOM absorption of PUR was computed to reduce 
f z GPP ez by as much as 24% at low solar elevations 
(75° S and day of year 260) when light was most limit- 
ing to phytoplankton growth (Fig. 4B). Reduction of 
f z GPP ez was less severe later in the year at higher sun 
elevations. This trend was not particularly sensitive 
to ozone concentrations ranging between 172 and 
344 DU (Fig. 4B, E). 

When the combined effect of CDOM absorption of 
PUR and UVR are integrated over the depth of the 
euphotic zone, it becomes apparent that the presence 


Fig. 3. (A) Seasonal changes in the mag- 
nitude of UV inhibition of gross primary 
productivity as a function of latitude for 
«3cix>m(450) = 0 and [0 3 ] = 344 Dob- 
son Units (DU). (B) Same as (A) but 
3cdom(450) = 0.03 m -1 . (C) Same as (A) 
but [0 3 | = 172 DU. (D) Same as (A) but 
[CDOM] = 0.03 nr 1 and [0 3 ] = 172 DU 



Day of the Year 
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Fig. 4. (A) Enhancement of the rate of 
depth-integrated daily gross primary 
production within the euphotic zone 
(j z GPP ez ) due to absorption of damag- 
ing ultraviolet radiation (UVR) by 
CDOM for [O 3 ] = 344 DU, calculated as 
the difference between corresponding 
simulations shown in Fig. 3A, B. (B) In- 
hibition of J Z GPP IV due to absorption 
of photosynthetically usable radiation 
(PUR) by CDOM for [0 :{ ] 344 DU. 
(C) Net inhibition of J^GPP^ due to the 
presence of CDOM at [O;*] = 344 DU. 
Includes both the enhancement and 
inhibition of J 7 GPP { . Z due to UV and 
PUR absorption, respectively, by 
CDOM. (D) Same as (A) but [ 0 :i ] = 172 
DU. (E) Same as (B) but [0 ;1 ] = 172 DU. 

(F) Same as (C) but [O;,] = 172 DU 



of CDOM results in a net decrease of J 7 GPP e7 (Fig. 4C). 
This is because the reduction in J z GPP oz due to the 
absorption of PUR by CDOM (Fig. 4B) was not fully 
compensated for by the slight photoprotection pro- 
vided by its absorption of damaging UVR (Fig. 4A). 
The net decrease in J z GPP e7 was ~19 to 24% during 
periods of low light levels in the presence CDOM 
(Fig. 4C), diminishing to 12 to 16% later in the year. 
Results were similar when ozone concentration was 
reduced to 172 DU (Fig. 4D, E, & F). 

Vertical profiles of downwelling irradiance and 
noontime production for 2 simulations representing the 
extremes in the irradiance field (August 28, 75° S and 
45° S, December 26) were examined to ascertain the 
range of variability produced by the presence of 
CDOM. Both irradiance (Fig. 5A) and production 
(Fig. 5C) decreased exponentially with depth at the 


Fig. 5. Simulated vertical profiles of noontime downwelling ir- 
radiance and production for August 28 at 75° S and for De- 
cember 26 at 45° S in waters with (w/) and without (w/o) 
CDOM. These simulations were chosen because they repre- 
sent the extremes in the range of solar zenith angles for our 
region of interest. The hatched areas labeled PP + and PP _ rep- 
resent the depth-integrated differences where production is 
enhanced and reduced, respectively 


log irradiance (pEin nv 2 s 1 ) Production (mg C rrr 3 h* 1 ) 
-24 -18 -12 -6 0 6 0 0.05 0.1 0.15 0.2 0,25 
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low light regime of 75° S on August 28. UVA, UVB, 
PUR, and production were reduced in waters contain- 
ing CDOM. This pattern was also evident at the high 
light environment of 45° S on December 26, although 
the absolute irradiances of all 3 wavelength bands 
were greater (Fig. 5B). Production, in this case, 
revealed a subsurface maximum (Fig. 5D) due to UV 
inhibition in surface waters. In the presence of CDOM, 
the depth of the production peak was shifted upward 
slightly (—10 m) relative to waters without CDOM, 
resulting in increased rates of production near the sur- 
face and reduced rates at greater depths (Fig. 5D). 

Calculating the total radiation absorbed by CDOM 
above a given depth illustrates that depth-dependent 
differences in production were due to changes in the 
vertical light field resulting from CDOM absorption 
(Fig. 6). Absorbed radiation was determined by inte- 
grating the difference in irradiance spectra calculated 
at a given depth with and without CDOM over UVR 
and PUR wavelengths. On December 26 (day of the 
year 360) at 45° S, at a depth of 10 m, this difference 
was 38.7 pEin m" 2 s' 1 for UVR and 121.6 pEin m~ 2 s 1 
for PUR (Fig. 6A). The ratio of UVRiPUR absorbed by 
CDOM was 0.307. At 75 m, this difference was reduced 
to 0.03 pEin m“ 2 s" 1 for UVR and 34.0 pEin nr 2 s 1 for 
PUR (Fig. 6C), decreasing the ratio of UVR:PUR ab- 
sorbed to 0.008. Therefore, the relative absorption of 
UVR versus PUR was >35 times greater at a depth of 
10 m than at 75 m. This same general trend was also 
observed for the low light simulation of August 28 at 
75° S (Fig. 6B, D). Consequently, at shallower depths, 
where CDOM absorbs proportionally more UVR than 


PUR, productivity will be enhanced by the presence of 
CDOM. At greater depths where UVR is scarce, the 
presence of CDOM will have little photoprotective 
effect. 

The depth-dependent pattern of primary production 
(PP) predicted in the presence and absence of CDOM 
(Fig. 5D) is useful for understanding how CDOM alters 
rates of productivity throughout the water column. The 
region of Fig. 5D labeled PP + represents those depths 
where the presence of CDOM enhances productivity. 
In this region of the water column, the UV photopro- 
tection provided by CDOM is greater than the reduc- 
tion in productivity due to the removal of PUR by 
CDOM. In contrast, PP" represents those depths where 
the presence of CDOM results in a net reduction in 
productivity. The size of the area of regions PP + and 
PP“ is a direct measure of the change in noontime 
depth-integrated production (mg C m -2 h" 1 ) due to the 
presence of CDOM. The percent change in primary 
production due to the presence of CDOM can be de- 
scribed by the quantity PP* (%), which is calculated as 

pp- = m.ioo <3) 

PP 

where PP is the rate of primary production when no 
CDOM is present. Values for PP* greater than 0 indi- 
cate a net increase in depth integrated production 
(j z GPP ez ) due to the presence of CDOM. Changes in PP* 
will be used below to evaluate the relative importance 
of various biotic and abiotic factors influencing J z GPP ez . 

Increasing a C DOM(450) always resulted in a lower 
PP* (Table 2). For example, an increase from 0.03 m" 1 
(as in the standard run) to 0.07 m" 1 resulted in a 39% 



Wavelength (nm) Wavelength (nm) 


Fig. 6. Subsurface 
spectral irradiance in 
waters with (w/) CDOM 
[acuoM(450) = 0.03 nr 1 ] 
and in waters without 
(w/o) CDOM at depths 
of (A) 10 m on Decem- 
ber 26, (B) 10 m on 
August 28, (C) 75 m on 
December 26, (D) and 
75 m on August 28 
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Table 2. Sensitivity of model to z riX)M (the depth over which CDOM is 
distributed) and <3 cdom( 450). PP + (mg Cm 2 h ’) denotes the noontime 
depth-integrated net increase in production at those depths where the 
presence of CDOM enhances productivity and PP (mg C m 2 h ] ) 
denotes the noontime depth-integrated net decrease in production at 
those depths where the presence of CDOM reduces productivity (see 
Fig. 5). PP* is the percent change in primary production due to the 
presence of CDOM [(PP + - PP }/PP, see 'Results']. z (:IX)M is the depth 
over which CDOM is distributed within the water column. Unless 
otherwise noted, variable values for all simulations were: day of year 
= 350, latitude = 45° S, [0 3 | = 344 DU, [chi a] = 0.1 mg m :i , 5 = 0.015 
nm 1 , and C:chl a = 50 (g:g) 


Parameter 


PP + 

pp - 

PP* 

<3cdom( 450) = 0.01 m 

\ z CD oM = 100 m 

0.35 

0.93 

-4.5 

3cdom( 450) = 0.03 m 

\ Zcdom = 100 m 

0.70 

2.60 

-14.8 

<3cdom( 450) = 0.05 m" 

\ Zcdom = 100 m 

0.89 

3.92 

-23.6 

tfcDOM(450) = 0.07 m 

\ Zcdom = 100 m 

0.97 

4.92 

-30.6 

a (T>o M (450) = 0.01 m 

\ Zcdom = 30 m 

0.35 

0.07 

2.2 

a< i X )m( 450) = 0.03 m 

\ zcdom = 30 m 

0.71 

0.78 

-0.6 

3( 3 )om( 450) = 0.05 m 

\ Zcdom = 30 m 

0.89 

1.82 

-7.2 

a ax)M (450) = 0.07 m 

\ Zcdom = 30 m 

0.97 

2.88 

-14.8 

a CIX)M (450} = 0.03 in 

\ Zcdom = 20 m 

0.74 

0.03 

5.5 

a CDOM(450) = 0.01 m 

\ Zcdom = 10 m 

0.78 

0 

6.0 

«3cdom( 450) = 0.03 m 

\ Zcdom = 10 m 

1.47 

0 

8.8 

<3cdom( 450) = 0.05 m 

\ Zcdom = 10 m 

1.57 

0 

10.1 

<3cdom( 450) = 0.07 m 

\ Zcdom = 10 m 

1.40 

0 

11.4 

3cdom( 450) = 0.03 nr 

\ Zf;DOM = 1 m 

1.72 

0 

13.4 


increase PP + and an 89% increase in PP . Consequently, 
PP* dropped from -14.8 to -30.6, indicating that the 
additional absorption by CDOM further suppressed 
rates of J z GPP oz . Furthermore, for all simulations where 
CDOM was distributed uniformly with depth, PP* was 
negative. This finding suggests that the decrease in 
J z GPP ez observed in the presence of CDOM in the stan- 
dard run is applicable over a wide range of CDOM 
concentrations. 

The magnitude of PP* was also a function of the ver- 
tical distribution of CDOM. As CDOM (with fixed ab- 
sorption properties, i.e. a CDO M(450) = 0.03 m" 1 and S = 
0.015 nm" 1 ) was increasingly restricted to the shal- 
lower waters, PP + increased while PP rapidly dropped 
to zero (Table 2). For example, when CDOM was 
restricted to the upper 30 m, PP + increased slightly (rel- 
ative to when CDOM was distributed over the 
euphotic depth) from 0.70 to 0.71 mg C nr 2 h" 1 while 
PP“ fell from 2.6 to 0.78 mg C m" 2 h _1 ; accordingly, PP* 
rose from -14.8 to -0.6. When CDOM was distributed 
only within the upper 20 m of the water column, PP + 
had increased to 0.74 mg C m 3 h 1 , PP had fallen to 
0.03 mg C m 2 h" 1 , and PP* increased to 5.5. PP + 
reached its peak of 1.72 mg C m“ 2 h" 1 when CDOM 
was restricted to the topmost 1 m of the water column. 
These results indicate that when CDOM is restricted to 
within 20 to 30 m of the surface, its presence will 
enhance j z GPP ez . 


When CDOM was distributed throughout the 
water column (the depth down to which CDOM 
is distributed within the water column will be 
referred to hereafter as z C dom), J z GPP e/ always 
declined. In addition, if z aX )M is <30 m, 
the presence of CDOM can enhance J Z GPP 0Z 
(Table 2). It follows then, that there exists a 
z C rx)M at which PP* =0, defined here as z ( | X)M *. 
When z ( IX)M > z rr)OM *, j z GPP ez declines in the 
presence of CDOM. Conversely, when z r DOM < 
Zcdom*- LgPP«z increases. 

The value for z CDOM * is controlled in part by 
the absorption properties of the CDOM; if the 
CDOM is strongly absorbing, then z c l)OM * must 
be shallower for CDOM to enhance J Z GPP, V . For 
example, when a rDOM (450) = 0.01 m 1 and 
Z(’dom wa s set in the model to 30 m, then PP* = 
2.2 (Table 2), indicating that the presence of 
CDOM enhanced j z GPP«, z and that z rDOM * was 
>30 m. (If z CDOM * were equal to 30 m, then PP* 
would have been 0; if Zctdom* were <30 m, then 
PP* would have been negative.) When 
acDOM(450) was increased to 0.03 m" 1 , then 
PP*= -0.6, demonstrating that z ri)0 M* was <30 
m. Within the range of values typically 
observed for a C DOM(450) in open waters, 
ZrDOM* is greater than 10 m. This can be seen 
for the case where z ax)M was set to 10 m (Table 2) and 
PP* was always positive. 

Increasing S, the slope of CDOM absorption spectra, 
increased PP*. For example, raising S from 0.012 to 
0.030 nm" 1 resulted in a rise in PP* from -17.2 to -5.0 
(Table 3). The increase in S caused a reduction in 
<3cdom(^) (see Eq. 2) at wavelengths longer than 
450 nm and an increase in a CDOM (?t) at shorter wave- 
lengths. Consequently, the change in PP* was due 
both to an increase in PP + resulting from the additional 
absorption of UVR by CDOM in surface waters and to 
a decrease in PP due to the reduction in the absorp- 
tion of PUR at depth. The depth-dependent change 
in production resulting from changes in S was small, 
however. 

Table 3. Sensitivity of model to the spectral slope, S, from 
Eq. (2). PP + (mg C m 2 h '), PP (mg C m 2 h ’), and PP* (%) 
are defined in Table 2. Unless otherwise noted, variable val- 
ues for all simulations were: day of year = 350, latitude = 45° S, 
Zcdom = 100 m - |0 3 ] = 344 DU, [chi a] = 0.1 mg m { , a riX)M (450) 
= 0.03 m \ and C:chl a = 50 (g:g) 


Parameter 

PP + 

pp - 

pp* 

S= 0.012 nm 1 

0.59 

2.80 

-17.2 

S = 0.018 nm 1 

0.82 

2.43 

-12.5 

S = 0.024 nm' 1 

1.05 

2.13 

-8.4 

S = 0.030 nm 1 

1.25 

1.89 

-5.0 
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Table 4. Sensitivity of model to chlorophyll a, PP* (mg Cm 2 
h 1 ), PP (mg C nr 2 h 1 ) F and PP* {%) are defined in Table 2. 
Unless otherwise noted, variable values for all simulations 
were: day of year = 350, latitude = 45° S, z riX)M = 100 m - [0 3 ] = 
344 DU, [chi a] = 0.1 mg m \ S = 0.015 nm \ a ( I)nM (450) = 
0.03 m \ and C:chl a = 50 (g:g) 


Parameter 

pp+ 

PP 

PP* 

[chi a] = 0.1 mg m 3 (standard run) 

0.70 

2.60 

-14.8 

[chi a] = 0.5 mg nr 3 

3.19 

11.57 

-14.4 

[chi a] = 1.0 mg m 3 

5.64 

19.05 

-13.1 

[chi a] = 5.0 mg nr 3 

14.39 

23.04 

-3.5 

[chi a] = 10.0 mg nr 3 

16.03 

15.70 

0.1 


Similarly, increasing concentrations of chi a through- 
out the water column will increase PP*. As chi a 
increased from 0. 1 to 5 mg m' 3 in the standard run, PP * 
increased from -14.8 to -3.5 (Table 4), indicating that 
the presence of CDOM will reduce depth-integrated 
primary production within a wide range of chi a con- 
centrations. Only at concentrations of chi a greater 
than -7 mg m" 3 did the addition of CDOM result in a 
net increase in depth-integrated primary production 
(PP* > 0) (Table 4). At high chi a concentrations, i.e. 
10 mg m 3 the addition of CDOM with an ^cdomM > 
0.015 m' 1 will drive PP* positive. 

The wavelength-dependent biological efficiency for 
damage to photosynthesis by UV, £(>.) [(pEin nr 2 s -1 )' 1 , 
Cullen et al. 1992] and the flux of UV radiation were 
used by the model to calculate UV inhibition of pri- 
mary production. Because relatively few measure- 
ments of e(A.) have been made, it was important to 
determine the sensitivity of the model to changes in 
e(X). Simulations showed that even when e(V) was 
increased 10-fold, PP* was only reduced by -65% and 
remained negative (-5.4) (Table 5). This indicates that 
if phytoplankton are as much as 10 times more sensi- 

Table 5. Sensitivity of model to e(k), the biological weighting 
function for the inhibition of phytoplankton photosynthesis by 
UV radiation (Cullen et al. 1992). PP + (mg C nr 2 h 1 ), PP (mg 
C m 1 h *}, and PP* (%) are defined in Table 2. Unless other- 
wise noted, variable values for all simulations were: day of 
year = 350, latitude = 45° S, z CUOM = 100 m, |0 3 ] = 344 DU, 
[chi a] = 0.1 mgm 3 S= 0.015 nm 1 , a ( -| X -,M{450) = 0.03 m" 1 , and 
C:chl a = 50 (g:g) 


Parameter 

PP + 

pp- 

pp . 

Standard run 

0.70 

2.60 

-14.8 

e(A) [(mW m' 2 )" 1 ] 

was increased by a factor of 2 

1.03 

2.58 

-12.7 

e{X) [(mW m 2 ) _1 ] 

was increased by a factor of 4 

1.38 

2.52 

-10.0 

e(jg [(mW m 2 ) x ] 

was increased by a factor of 10 

1.89 

2.43 

-5.4 


tive to the presence of UVR than currently available 
photophysiological data suggest (e is high), the addi- 
tion of CDOM will still reduce depth-integrated pri- 
mary production. 

DISCUSSION 

The relationship between the presence of CDOM 
and J Z GPP«, Z is a complex one, depending upon the 
absorption properties and vertical distribution of 
CDOM (Table 6). A few generalizations are possible, 
however, based upon model results. Because UVR is 
most abundant within a few tens of meters of the sea 
surface, it is at these depths that light absorption by 
CDOM has the greatest net positive influence on rates 
of primary production, acting as a UV photoprotectant. 
In deeper waters, where little UVR remains, the pres- 
ence of CDOM can only reduce rates of production via 
the absorption of PUR necessary for carbon fixation. 
Therefore, the balance between absorption of UVR 
and PUR by CDOM will determine whether the pres- 
ence of CDOM enhances or inhibits j z GPP ez . 

When both CDOM and phytoplankton are distrib- 
uted uniformly within the euphotic zone, the presence 
of CDOM always reduces J z GPP ez , regardless of the 
value of <3 c;dom( 450). The only exception appears to be 
when chi a concentrations exceed ~7 mg m 3 (Table 4), 
a condition rarely observed in the open ocean. In 
waters with more typical concentrations of chi a, the 
amount of UV photoprotection provided to phyto- 
plankton by CDOM in surface waters is not as great as 
the reduction in the depth integrated rate of primary 
production due to the absorption of PUR. This pattern 
was found to be independent of either the magnitude 
and spectral shape of ScdomM or the sensitivity of the 
phytoplankton to UVR. For example, when e(^) was 
increased 2- to 10-fold (phytoplankton were rendered 
highly sensitive to UVR), the presence of CDOM still 
reduced j 2 GPP ez , suggesting that the predicted de- 
crease in j z GPP ez was not merely the result of an unre- 
alistically low biological weighting function. 

Table 6. The effects of z CDOM and a CIXJM (450) on integrated 
gross primary productivity (j z GPP ez ). The following gen- 
eralizations apply only when [chi a] <7 mg nr 3 and 
S- 0.015 nm 1 


z CDOM ( m ) 

acDOM(450) (m ] ) 

L GPPez 

z CDOM > 30 

<0.07 

Reduced 

10 < z CD om <30 

>0.015 J 

<0.015 ] 

1 

r Enhanced 

Z CDOM —16 

<0.07 J 

I 
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Unfortunately, relatively little field data are avail- 
able describing the vertical distribution of CDOM 
within the water column. Our results suggest that 
regions characterized by deep mixed layers and 
approximately uniform chi a profiles, such as the North 
Atlantic or the open Southern Ocean, should experi- 
ence a decline in J z GPP ez when CDOM is present. The 
Southern Ocean is of particular interest in this regard 
because of its exposure to high fluxes of UVB during 
periods of reduced stratospheric ozone concentration. 
Numerous field studies have shown that phytoplank- 
ton in surface waters experience reduced rates of 
photosynthesis when exposed to naturally occurring 
levels UVB (Karentz et al. 1991, Helbling et al. 1992, 
Smith et al. 1992, Holm-Hansen et al. 1993a, b, Cullen 
& Neale 1994). Unfortunately, the influence of CDOM 
was never investigated during these studies to deter- 
mine whether its presence was influencing the mea- 
sured rates of UVR photoinhibition. Our model pre- 
dicted that while under conditions of increased UVB 
flux (i.e. when ozone concentrations were reduced 
by 50 %) photoinhibition in near-surface waters was 
reduced when CDOM was present, but | z GPP oz contin- 
ued to decline, even at high latitudes. In short, the 
presence of CDOM in the open Southern Ocean would 
probably diminish, but not eliminate, decreases in 
j y GPP ez associated with increased UVB fluxes beneath 
the ozone hole. 

The presence of CDOM should enhance j z GPP ( . z in 
areas of equatorial upwelling, marginal ice edge 
zones, frontal zones, and eddies. These are locations 
where <3 C dom(450) is likely to be small (less than 
-0.015 m 1 ) or z CDOM shallow (less than -40 m) (Table 2) 
or where phytoplankton are restricted to near-surface 
waters. J^GPP^ will be enhanced under these condi- 
tions because even CDOM with a relatively low ab- 
sorption coefficient will absorb a substantial fraction of 
UVR at the sea surface relative to its absorption of 
PUR. These conditions favor enhanced phytoplankton 
growth. As CDOM absorption increases, however, the 
benefits of UVR absorption will diminish because pro- 
portionally less UVR will be left to be absorbed at 
a given depth and the availability of PUR will be 
reduced, causing a reduction rather than an increase 
in j z GPP ez . Similarly, when z c:DOM increases, little UVR 
will be transmitted to depth (it is effectively absorbed 
in surface waters by CDOM and seawater), and the 
presence of CDOM will reduce PUR in deeper waters 
and / z GPP ez will decline. 

The beneficial effect of CDOM absorption on J z GPP ez 
will be exaggerated when the spectral slope, S, is large 
because the ratio of absorbed UVR:PUR will increase. 
This may be an important consideration in the open 
ocean where values for S tend to be higher (0.015 to 
0.029 nm" 1 ) than those measured from CDOM col- 


lected in coastal regions (Bricaud et al. 1981, Green & 
Blough 1994). Therefore, the presence of CDOM is 
most advantageous to depth-integrated phytoplankton 
production when a rDO M(450) is small, S is large, and 
z ct>om < 30 m. 

When CDOM is highly absorbing, the effect of CDOM 
on J z GPP ez depends upon a CDO M(450), z n)OM , and the 
depth over which phytoplankton are distributed (z rhl ri ). 
This is the most difficult case to generalize. Clearly, 
however, as a ( I)O m(450) increases, z rI)OM or z rh id must 
decrease accordingly (to -30 m) in order for the pres- 
ence of CDOM to enhance j z GPP e/ . CDOM and phyto- 
plankton are often restricted to the surface layer in 
waters that are markedly stratified and exhibit a high 
degree of biological activity, such as in coastal regions 
during the spring and at river mouths. CDOM concen- 
trations >2 g CDOM-C m :i have been reported in 
coastal regions such as the Delaware Bight, USA, 
(DeGrandpre et al. unpubl.) and in the Orinoco River, 
Venezuela, outflow plume (-3.1 g CDOM-C m *; 
Blough et al. 1993). At concentrations of this magni- 
tude, even weakly absorbing CDOM would tend to 
remove a great deal of the PUR available to phyto- 
plankton and cause a reduction in J z GPP t , z . However, 
because terrigeneous sources of CDOM generally have 
a high specific absorption coefficient [a Y DOM (450) = 
0.51 to 1.76 m 2 g CDOM-C m Zepp & Schlotzhauer 
1981, Blough et al. 1993, and Blough & Green 1995), 
phytoplankton associated with river outflow plumes 
would experience reduced depth integrated rates of 
production unless either z rDO M or z rh | rf was very small 
(<10 m ). 

Reducing [0 3 ] in the model had little influence on 
the magnitude of the change in J z GPP t?z resulting from 
light absorption by CDOM. This insensitivity to [0 3 ] 
was predicted regardless of a CT)O M(450) and whether 
CDOM was distributed uniformly throughout the 
euphotic zone or restricted to surface waters. This was 
somewhat surprising considering that the flux of UVB 
at the sea surface is much greater when [0 3 ] is low. In 
fact, under conditions of reduced ozone the increased 
flux of UVB at the sea surface reduced PP + slightly. At 
depth, PP - remained relatively unchanged. Conse- 
quently, while PP* tended to decrease when [0 3 ] de- 
creased, it did so only slightly. Therefore, even in 
waters where the flux of UVB is relatively high, such as 
the Southern Ocean when the ozone hole is present, 
CDOM is expected to provide little or no added bene- 
fit to phytoplankton production. 

The effect of differences in phytoplankton commu- 
nity structure on j z GPP 0/ was not considered but may 
have important implications when applying model 
results to field conditions. For example, Karentz et al. 
(1991) have shown that interspecific differences in 
UVB sensitivity among Antarctic diatoms can vary by a 
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factor of <50. However, the value for e(>.) determined 
for a single species, Phaeodactylum sp., was used in 
the model to estimate UV inhibition. If this diatom spe- 
cies proves to be particularly sensitive to UVB (relative 
to other species), then the presence of CDOM will be 
less beneficial than predicted by the model, even in 
surface waters where the flux of UVB is greatest. The 
use of incorrect values for e(A.) is not likely to change 
our general conclusions, however, because the model 
was not particularly sensitive to e(A,) (Table 5). 

Although the optical properties of CDOM have 
received a great deal of attention in recent years, little 
is currently known about vertical distributions of 
CDOM in the water column and how absorption prop- 
erties vary with depth or with season. Even less is 
known about the dynamics of CDOM degradation and 
formation in the open ocean. Although remote sensing 
techniques have been used recently to measure 
CDOM absorption over wide areas of the ocean (Hoge 
et al. 1995a, b, Vodacek et al. 1995), much more must 
be learned about the distribution of CDOM on a global 
scale before the results of the present study can be 
generalized to a wide range of field conditions. This is 
particularly true for regions of high biological activity 
where the precise effect of CDOM on j z GPP ez is deter- 
mined by z C[ )om and a C DOM(450), quantities which are 
virtually unknown for most waters, and by z chld . Even 
in deeply mixed water columns where the presence of 
CDOM always reduces J z GPP ez , the magnitude of this 
reduction cannot be estimated without knowledge of 
CDOM distributions and optical properties. 
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